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// is proved that the integral equation 

f G(x)F(xy)H(y)f(y) dy = \f(x) 

has at least one nonzero eigenvalue if F is any integral junction of finite 
order, G and H are any bounded functions on [ — 1,1], and the trace of the 
kernel G(x)F(xy)H(y) does not vanish. In particular, this theorem furnishes 
the first rigorous proof that the kernel exp [ik(x — y) 2 ], which arises in the 
theory of the gas laser, has an eigenvalue for arbitrary complex k. 

I. INTRODUCTION AND SUMMARY 

In an idealized model of the gas laser or optical maser, as studied by 
Fox and Li 1 - 2 and others, electromagnetic radiation is reflected back and 
forth between two infinitely long metal strips which are mirror images of 
each other. A typical field quantity, such as the current density, at the 
surface of each reflector satisfies the integral equation 

jf exp{i[k(x - yf - h(x) - h(y)]\ f(y) dy = \f(x), (1) 

where k is a dimensionless real parameter which depends on the width 
and spacing of the reflectors and the wavelength, and h(x) is a real func- 
tion specifying the departure of the reflecting surfaces from parallel 
planes. 

The eigenf unctions of (1) represent the field distributions at the re- 
flectors of the possible modes of oscillation of the laser, and the eigen- 
value X corresponding to a particular mode represents the complex factor 
by which the field strength is multiplied as a result Of one reflection and 
transit between the reflectors. From the magnitude of X one can deduce 
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the amount of amplification which would have to be provided by an ac- 
tive medium between the reflectors in order just to sustain oscillations 
in the given mode, while the phase of A determines admissible reflector 
spacings for oscillations at a particular frequency. 

The mathematical interest of (1) centers around the fact that its 
kernel K(x,y) is complex symmetric but not Hermitian;* that is, 

K(x,y) = K(y,x) but K(x,y) * K(yjj- < 2 ) 

The ordinary theory of Hermitian kernels does not even suffice to prove 
the existence of eigenvalues of complex symmetric kernels. Fox and Li 1 
have made extensive calculations of the eigenvalues and eigenfunctions 
of (1) for h(x) = by iterative numerical techniques up to about k = 60 
(in applications k may be as large as a few hundred) ; but heretofore there 
has been no formal mathematical proof of the existence of solutions 
except f for | k | « 1, which is not a case of physical interest. 

This paper contains a proof of the following 

Theorem: Let G(x) and H(x) be any bounded functions on the interval 
— 1 ^ X ?=■ 1, and let F{z) be any integral function of finite order such that 

f G(x)F(x 2 )H(x) dx * 0. (3) 

Then the integral equation 

G(x)F(xy)H(y)f(y) dy = \f(x) (4) 



/ 



has at least one nonzero eigenvalue. 

As a corollary, it follows that the integral equation (1) has at least one 
eigenvalue for arbitrary complex k, provided only that 



/: 



£f 2 ' A(x) dx 9* 0. (5) 



Furthermore if h(x) is an even function of x, then (1) has at least two 
eigenvalues for all but certain exceptional values of k, a particular excep- 
tional value being k = 0. 

The idea of the proof is quite simple. The assumption that F{xy) in 
(4) is an integral function of finite order means that ultimately the coeffi- 
cients of its Taylor series in powers of xy fall off with extreme rapidity. 

* The kernel is normal in the special case h(x) = kx 2 . The eigenfunctions of 
exp (— 2ikxy) are prolate spheroidal wave functions, as pointed out in connection 
with lasers by Boyd and Gordon. 3 

t If | k J « 1 then exp [ik(x - r/) 2 ] is nearly unity, and the existence of at least 
one eigenvalue follows from perturbation theory; see Sz.-Nagy.< 
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If we truncate the Taylor series after a finite number of terms, (4) is 
replaced by an integral equation with a kernel of finite rank. The eigen- 
values of such a kernel are merely the latent roots of a finite matrix, and 
these are not all zero if their sum, which is the trace of the matrix, does 
not vanish. The limiting value of the trace is just the left side of (3), 
and does not vanish by hypothesis. By taking more and more terms of 
the series for F<xy), we obtain a sequence of larger and larger matrices, 
whose elements ultimately vanish very rapidly with distance from the 
upper left corner. We show that it is possible to pick one eigenvalue from 
the set of eigenvalues of each succeeding matrix in such a way that the 
resulting sequence of numbers has a nonzero limit point. This limit point 
is an eigenvalue of the infinite matrix, and hence an eigenvalue of the orig- 
inal integral equation. 

Details of the argument just sketched are given in a series of lemmas 
in the next section, followed by the proof of the main theorem. Since the 
existence proof makes heavy use of asymptotic inequalities, it does not 
generally provide a practical technique for obtaining numerical results. 
The important practical question of finding approximate expressions, 
valid for large A - , for the eigenfunctions and eigenvalues of equations such 
as (1) is a separate problem, as is also the question whether any particular 
equation has a finite or infinite number of eigenvalues. 

For a gas laser with finite (not strip) mirrors of arbitrary, dissimilar 
shape and size, the integral equation still has a complex symmetric 
kernel, 2 although the domain of integration is two-dimensional and the 
kernel is more complicated than that of (1). The existence of eigenvalues 
in the most general case still remains to be settled. 

II. MATHEMATICAL DETAILS 

We shall use the following notation referring to an n X n matrix: 
A M = (atj), i= 1,2, ■■■ ,n; j = 1, 2, ■ • ■ , n; 

-1""(0 = Z I «■•;!, * = 1,2, ■■■ ,n; (6) 

S(A (n) ) = ±A^(i) = ZZ|«,,|. 
i=i i-i j=i 

If the superscript is omitted, n is understood to be infinite. 
Lemma 1: 



detA l " } \ ^ UA (n) (i). (7) 
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Proof: Using Hadamard's inequality, 



| det A {n) | ^ II 



E I a « 



=i _y=i 



1/2 



sn 

t=i 
Lemma 2: 
det(A M + B M ) -detA {n) 



E I <m 



2H1/2 n 



(8) 



^ n^ w w + & n) (i)] - n^ (n) (0- 



(9) 



i=i 



Proof: The lemma is obviously true for n = 1. To proceed by induction, 
assume it is true for all determinants of order n — 1, and expand the 
determinants in (9) by minors of the first row. Let &,- be the algebraic 
complement of thj + hi in A (n) + B M , and let 4i, be the algebraic com- 
plement of aiy in A {n) . Then 



det(A M + fl ( ") = E («w + &ii)Ci> 



./ -i 



(10) 



- det A (n) + E ffly(C w - At/) + E bifiu . 



By Lemma 1, 



(11) 



! C„ | ^ fl [E I a« + 6* |] 
t=2 |_*=i J 

i=2 

By the inductive hypothesis, 

\C U - Au\ £ fl[A M {i) + B ln \i)) - n^ (B) (»). (12) 

where we have used the fact that the right-hand side is increasing as a 
function of the A (n) (i) and B {n \i). Hence (10) gives 



EXISTENCE OF EIGENVALUES 117 

\det(A (n) + B {n) ) - detA (n) | 

^ A (n \l) lfl[A <n \i)+B M (i)] - U.A M (i) 

[i=2 i=2 

+ B'"\l)fl[A< n) (i) + B M (i)] (13) 

i=2 

= Il[A M (i)+B ln) (i)]-ILA M (i), 

and the induction is complete. 

Now let (B be the Banach space* whose elements are all bounded se- 
quences of complex numbers, e.g., 

x = (xi , x 2 , ■• ■ ,x iy •••) (14) 

with norm 

|| x || = sup \xi\. (15) 

i 

Let A be a linear matrix operator on the space ($>, defined by 

(Ax)i= f^aijxj, i= 1,2, ••-. (16) 

Ax will be an element of (B provided that sup A(i) is finite. The norm of 
A is defined by 

|| A || = sup { || Ax ||; || a; || = 1}, (17) 

and it is easy to show that 

|| A || = sup^l(t). (18) 

t 

Henceforth we shall restrict our attention to matrix operators for 
which 

S(A) m 'EA(i) < oo. (19) 

1=1 

Such operators are completely continuous, because they can be ap- 
proximated by the sequence {A " } of completely continuous operators 
which converges in norm to .4. Here A (n) is a matrix whose elements co- 

* The standard definitions and theorems which we shall require from functional 
analysis may be found in Kolmogorov and Fomin. 6 
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incide with those of A for 1 ^ i ^ n and 1 ^ j ^ n, and are zero other- 
wise. 

A complex number X is said to be in the spectrum of an operator A if 
the operator A — XI has no inverse. An eigenvalue of A is any value of 
X for which there exists a nonzero x satisfying the homogeneous equation 

Ax - Xx = 0. (20) 

If A is completely continuous and if X ( ^ 0) lies in the spectrum of A, 
then X is an eigenvalue of A, In finite-dimensional space the eigenvalues 
are the latent roots of the matrix A (B) ; that is, they are the roots of the 
characteristic equation 

det (A 00 - X/ (n) ) = 0. (21) 

Lemma 3: If A {n) has X as an eigenvalue, then A + B " has X', 
where 



X-X'\ ^ lf[[A ln \i) + B M (i) + |X|] 



-n^ ( " } (0 +1*11 



l/n 



(22) 



(n) ■ D («) 



Proof: Denote the eigenvalues of A " + B " by Xi , X 2 , • • • , X„ . 
Then 



(X - Xi)(X - X 2 ) ••• (X - x„) I 

= I det (A in) + B (n) - Xl (n) ) - det (A {n) - X/ ( "> 



(23) 



the second determinant being equal to zero because X is an eigenvalue of 
A {n) . Let 

so that 

D (n) (i) = t, I <*a - W« I ^ ^ <B) (*) + I X |. (25) 

Then, using Lemma 2, 

n i x - x* i ^ ri p w (*) + ^ (n) (^] - n D (n \i) 
k=i i=i «=i 

^nU (n, W + * (B) (0 + |X|] (26) 

i=l 

-n^ (n) (*) + ix|], 
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since the right side of the first line is an increasing function of D " (i). 
It follows from (26) that for at least one of the factors | X — X* | the 
inequality (22) holds. 

Lemma 4- Let A be an infinite matrix with S(A) < °o. Suppose that 
from the eigenvalues of the sequence of finite matrices [A ln) ] we can pick a 
sequence {X (n) } such that X*' does not approach zero as n — > °o . Then A has 
a nonzero eigenvalue. 

Proof: The X (n) are bounded, since in fact 

| X (n) | ^ || A (n) || = max A (n \i) ^ S(A). (27) 

Also for sufficiently large n we can pick a subsequence which is bounded 
away from zero, and which therefore has at least one nonzero limit point. 
Suppose that the subsequence X converges to the limit point X t^ 0, 
as p runs through some increasing sequence of integers. We assert that 
X is an eigenvalue of A. If it were not so, then (A — \I)~ would exist 
and therefore be bounded. Suppose (A — \I)~ l were bounded, and let 
x (p) be the characteristic vector of A p corresponding to X p . Then we 
would have 

x ™ = (A - \l)~\A - *I)x (p) 

= (.4 - x/r 1 [,i (p V p) - x (p V p) 

(28) 
+ (A - A (p) )x ip) - (X - \ (p) )x (p) ] 

= (A - XI)~ : [(^ - A (p) )x <p) - (X - X (p) )x (p) ], 

where in the last equation A {p) represents an infinite matrix which coin- 
cides with A in a square of side p in the upper left corner, and has zeros 
elsewhere. Taking norms, we have 

|| x' p) || ^ \\(A - xir 1 !! \\(A - A (P V P) - (x - x (p V p) II 

^ || (i4 - XT") -1 || [|| A - A (p) || + | X - X (p) |] || x (p) ||, (29) 

or 

II U - X/)-' II i ii^^V+IX-A-l • ^ 

But since both || A — A (p) || and | X — X (p) | go to zero as p — > °°, we 
derive a contradiction. 

Theorem: Let A be an infinite matrix with S(A) < °o and with Tr(A) 
*0.If 

S(A) - S(A 0,) ) < (c/nT, (31) 

for some c, e > 0, then A has a nonzero eigenvalue. 
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Proof: Since Tr(A) ^ and Tr(A (n) ) -> Tr(A), it follows that for 
n ^ n\ (say) and some 5 > 0, we have | Tr(A (n) )| ^ S. Since the trace 
is the sum of the eigenvalues, A (n) must have at least one eigenvalue 
X (n) such that 

I X (n) | ^ S/n. (32) 

We shall in fact show that if n x is a sufficiently large fixed integer, and if 

iii-^m, .7 = 1,2,3, ■•• (33) 

then for each j there exists an eigenvalue which is uniformly bounded 
away from zero, i.e.. 

X ( "' } ^ S/2m . (34) 

Then by Lemma 4 the theorem will be proved. 
We substitute into Lemma 3 as follows : 



n = n j+1 , 

(35) 



\ I = I X ( "'° \ = t, 



A M = AW , 

b m = A (n *+ l) - A 1 ** , 

where it is understood that A (nj) now represents the original matrix 
A lnf) augmented below and to the right with enough zeros to give it 
dimensions fty+i X »y+i ■ Then (22) becomes 



g ("ft 1 [lM(i) + t ] - f'^- n 'f[[A (n '\i) + t] 



l/n y+ , 



(36) 



Since 

| x (»/> _ x (»*+i) j ^ f - | X (n ' +l) |, (37 

we can rearrange (36) to get 

J TT 



X («*+») | ^ f _ m [A(i) + 1] - r >+'-"'• nu (n '°(*') + <] 

I. i=l t=l 



(38) 



Hence 
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IX 



^*Hff[ i+ ^]-fi[ i+ ^]T 



*i- n 



i + 



y A <^ ' > (^)^\ ,/n '' +, 



^]-ft[ 1 + =«i>]} 



(39) 



since we already know that / ^ d/n } . 
Now consider 



nhfl^B^f],^ 



su) 



(40) 



Also 






»=„, + l 







(41) 



provided that th and hence ?i> are sufficiently large, where in the next 
to last step we have used (31) and in the last step we have used e x 
^ 1 + 2.x- for | a:^ 1, say. Finally, 

ft [ 1+ M^)] 



»fl[ 



1 + 



•ffl 



1 + 



n,vl(z) 



A = l 



»iA(i] 



1 + 



njA(i)' 



(42) 



> 



- '\ \su) - .s(,i ( "^)i fl [i + *££>] 
'-ffe)'"}n[-^]- 
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Substituting (40), (41), and (42) into (39) yields 

I X""+" | > . _ K(2 + B| ) ( ± y _ nsSU)Y ln " 



= 1 - 



V n,(2 + % ) -p 



(43) 



n," 2 



where in the last step we have used the fact that ra,-+i = 2wj and have 

set 

Cl = c l/2 exp[,S(A)/2*]. 
If we assume in advance that 

5^2, ni ^ max (2, 4/e), 
then 



L 5 J n/' 2 ~ L 5 J n/' 2 



2ci 



where 



c 2 = 



2ci 



Hence (43) and (46) imply 



~ &nfl* 2 { i-» t ' 4 



r = 2~ lli < 1. 



C2 « J- 1 

= c 2 r , 






A w^ 1 - «■"• 



>-ii 



(44) 

(45) 



(46) 



(47) 



(48) 



(49) 



and by induction 

I . (nj) | J-\ 

^r-j ^ n [i - *r - 

I X (ni) I >=^r 
But if c 2 ^ M, say, then 

ft (1 - c 2 r i_1 ) = exp [fj log (1 - c 2 r J_1 )J 

^ exp -2 2 c 2 r J_1 = exp I -j-~ 

where the last step requires 

< y 2 {\ - r) log 2 = y 2 {\- 2~ l,i ) log 2, (51) 



(50) 



>H, 



ca 
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and by (47) this inequality can always be satisfied for large enough n\ . 
But (49) and (50) imply 

X (B '° ^ V 2 \ (ni) ^ «/(2ni) > (52) 

for all j, and so the theorem follows from Lemma 4. Q.E.D. 

An integral function of finite order p is a function F(z) which has no 
singularities in any finite region of the z-plane, and whose maximum 
modulus M(r) on the circle \z\ = r satisfies 

log M (r) < r k (53) 

for all sufficiently large r when A; > p, but not when k < p. Such a func- 
tion may be expanded in a Taylor series, 

F(z) = Y, a„z", (54) 

n=0 

which converges for all z, and whose coefficients satisfy 6 

I On | < l/n nt (55) 

for all sufficiently large n, where e is any fixed number less than 1/p. 
Alternatively, for any fixed e < 1/p, there exists a constant c such that 
for all n > 

We are now ready to prove the result stated in Section I. 

Theorem: Let G(x) and H(x) be any bounded functions on the interval 
— 1 ^ x ^1, and let F(z) be any integral function of finite order such that 

J_ G{x)F{x 2 )H(x)dx ^ 0. (57) 

Then the integral equation 

f G(x)F{xy)H{y)f(y)dy = \f(x) (58) 

j- 1 

has at least one nonzero eigenvalue. 

Proof: Expand F{xy) in a Taylor series, so that the integral equation 
becomes 

/ Z [an-i m G(x)x'-% l ^ m H(y)y''- 1 ]f(y)dy = */(*). (59) 
Let 
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00 

fix) = G(x) E fnOn-lX n -\ (60) 

n=l 

where {/„} is a bounded sequence of complex numbers; the a„'s tend to 
zero fast enough so that f(z)/G(z) will be an integral fimction of finite 
order. 

Since the powers of x are linearly independent, (59) is equivalent to 
the matrix equation 

4f = V, ( 61 ) 

where 

aij = aji = (a.-xa,-!) 1 ' 2 f G{t)H{t)t i+i ~ 2 dt, (62) 

i = 1, 2, • • • ; j - 1, 2, • • • . 

Since G(x) and ff(.r) are bounded in -1 ^ a; ^ 1 and the Taylor co- 
efficients of F(z) satisfy (56), it is clear that 

In preparation for an application of the preceding theorem, consider 

= »!, [(I)"' Sr^ (?)"} 

Now (c/f) 3 ' 2 is bounded as j -> », and 

t_l_- a f-^ = lo g r 2 f^l, (65) 

pi « + j - 1 ■'i-i a; _ % - 1 J 

which is bounded for t ^ w + 1 ^ 2. Hence with a new bounding con- 
stant we have 

8(A) - S(A 0,) ) ^ Mx E (£X (66) 

*=n+l V / 

Choose log n ^ (2 + log c)/e, so that w € ^ ce 2 ; then 

oo / l/2\i /-oo / l/2\z /•«> / l/2\z 

.5. W*/.fe) **/.&■)* 

(c/n ) 



(log c — € log n)/2 



- < C^V 

2 = Vn" 2 / ' 



(67) 
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and so from (66) 

S(A) -S(A {n) ) ^fe)". (68) 

where Ci is a new bounding constant and ei = e/2. 
Finally we have 

Tr U) = Z a ti = E Oi-i f G(t)H(t)t 2i - 
= \\j{t)H{t)F(i 2 )dl, 



dt 



(09) 



and this does not vanish by hypothesis. Hence all the conditions of the 
previous theorem are satisfied, and the integral equation has a nonzero 
eigenvalue. Q.E.D. 

Since exp ( — 2ikz) is an integral function of finite order 1, it is an obvi- 
ous corollary that the kernel exp i[k(x— y) 2 — h(x)— h(y)] has a non- 
zero eigenvalue for arbitrary complex k, provided only that h(x) is 
bounded and that 



i: 



-2ih(x) 



e~ MW dx * 0. (70) 

Furthermore if h(x) is an even function of x and if f(x) is an even func- 
tion which satisfies 

f exp [i\k(x* + y 1 ) - h(x) - h(y)}\ cos (2kxy)f(y)dy = V 2 Xf(x), (71) 
•'o 

then f{x) also satisfies (1). But the theorem just proved obviously holds 
for arbitrary finite limits of integration and applies to the kernel of (71), 
so (71) has at least one nonzero eigenvalue if 

/ exp \2i\kx 2 - h(x)]} cos {2kx 2 )dx ■£ 0. (72) 

Similarly if h{x) is even and if f(x) is an odd function which satisfies 

I exp {i[k(x 2 + if) - h(x) - h(y)]} sin (2kxy)f(y)dy = JflX/fcr), (73) 

then f(x) also satisfies (1), and (73) has at least one nonzero eigenvalue 
if 



f 



exp [2i[kx 2 - h\x)}\ sin (2kx 2 )dx ^ 0. (74) 
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At least one of (72) and (74) will be satisfied whenever (70) holds. Except 
for certain particular values of k, one of which is evidently k = 0, both 
(72) and (74) will be satisfied, and (1) will have at least two distinct 
eigenf unctions corresponding to nonzero eigenvalues. 
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